SHARP SPECTRAL BOUNDS ON STARLIKE DOMAINS 



R. S. LAUGESEN AND B. A. SIUDEJA 

ABSTRACT. We prove sharp bounds on eigenvalues of the Laplacian that complement the 
Faber-Rrahn and Luttinger inequalities. In particular, we prove that the ball maximizes the 
first eigenvalue and minimizes the spectral zeta function and heat trace. The normaliza- 
tion on the domain incorporates volume and a computable geometric factor that measures 
the deviation of the domain from roundness, in terms of moment of inertia and a support 
functional introduced by Polya and Szego. 

Additional functionals handled by our method include finite sums and products of eigen- 
values. The results hold on convex and starlike domains, and for Dirichlet, Neumann or 
Robin boundary conditions. 



1. Introduction 

How do eigenvalues of the Laplacian depend on the shape of the domain? We will obtain 
new quantitative estimates on the eigenvalues in terms of explicitly computable geometric 
functionals. 

Write Xj for the Dirichlet eigenvalues of the Laplacian on the bounded domain f2 in 
M. d , d > 2, with corresponding L 2 -orthonormal eigenfunctions Uj, so that 

— Auj = XjUj in Q 
uj = on dil 

and 

< Ai < A 2 < A 3 < . . . . 

These eigenvalues represent physical quantities such as frequencies of vibration, rates of 
decay to equilibrium in diffusion, and energy levels of quantum particles. 

The problem of understanding how eigenvalues are affected by the shape of the domain, 
and of identifying domains that extremize eigenvalues, is long-standing and difficult. Nu- 
merous monographs and survey articles summarize the state of knowledge in euclidean 
space [3, 4, 7, 19, 23, 24]. Important results have been obtained on closed surfaces too, for 
example, see [22, 39]. 

Let us first describe our main result in the special case of 2 dimensions (Theorem 1.1). 
Later we extend to all dimensions (Theorem 3.1) and to Neumann and Robin analogues 
(Theorem 3.3 and Theorem 3.5), and finally to the sloshing eigenvalues (Section 13). 
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Figure 1 . A starlike domain with outer normal N(x). 



Consider a starlike domain as in Figure 1 and define scale-invariant geometric factors 

Go = f / -4^^)' G 1= 2 -^, (1) 
2vr J m x ■ N(x) A 2 

where N(x) is the outward unit normal vector, A is the area of fi, and /origin = J n \x\ 2 dA 
is the polar moment of inertia about the origin. Note that x ■ N(x) > because the domain 
is starlike. In higher dimensions we will later define Go and G\ differently, although the 
definitions will reduce to (1) in 2 dimensions. 

Let G = max{Go, Gi}. Then G > 1 for all starlike domains with equality for centered 
disks, by Lemma 2.2 below. Thus one may regard the value of G as measuring the deviation 
of the domain from roundness. This deviation can arise in two ways: a highly oscillatory 
starlike boundary would make G large, while an elongated boundary (such as an eccentric 
ellipse) would make G\ large. 

Now we can state the main result in the plane, where the disk maximizes eigenvalues of 
the Laplacian. 

Theorem 1.1 (Dirichlet in 2 dimensions). Suppose the function R{6) is 2ix -periodic, posi- 
tive, and Lipschitz continuous, and consider the starlike plane domain Q = {re %9 : < r < 
R{0)}. Let $ : K + — > M be concave and increasing. Then X]j=i ®(^jA/G) is maximal 
when Q is a disk centered at the origin, for each n > 1. 

Hence each of the following eigenvalue functionals achieves its maximum value when 
the domain Q is a centered disk: 

XxA/Go, X2A/G0, (\t + --- + X s n ) 1/s A/G, ^/X 1 X 2 ---X n A/G, 

for each exponent < s < 1. 

Further, the partial sums of the spectral zeta function and trace of the heat kernel are 
minimal when Q is a centered disk. That is, the functionals 

n n 

^2(XjA/G) s and ^ exp(-A i vit/G) 
i=i i=i 
attain their smallest value when Q is a centered disk, for each s < < t. 

It is natural in this theorem to multiply Xj by A, because Xj scales like 1/A. (Intuitively, 
low frequencies come from large drums.) 



3 



rotation U 




Figure 2. A linear- on-rays transformation from a domain of area tx to 
the unit disk. To insure that the mapping preserves area locally, we require 

RiQf de = d<f). 



The result for A 2 follows immediately from the one for Ai, because the ratio A 2 /Ai is 
maximal on the disk by Ashbaugh and Benguria's "sharp PPW inequality" [2]. 

Note the theorem is better for the first and second eigenvalues than for other functionals, 
in the sense that we normalize \iA/G and \ 2 A/G with G instead of G, and G < G by 
definition. 

The theorem improves on the standard "inradius bounds" for Ai and \ 2 , on convex do- 
mains. To see this, write D in for the largest open disk centered at the origin and contained 
in f2. Then Aj(Sl) < Aj(D in ) for all j, by domain monotonicity of the Dirichlet spectrum. 
Theorem 1 . 1 implies this inradius bound for j = 1 , 2, as one sees by using that A in < A/Go 
(Lemma 10.2); here A in is the area of D in . 

Our theorem significantly extends the only known result of its type, which is the case 
(n = 1) of the fundamental tone Ai with Dirichlet boundary conditions. That case is due 
to Polya and Szego in 2 dimensions and Freitas and Krejciffk in higher dimensions, as 
explained after Theorem 3.1. 

To treat higher eigenvalues, we need a fundamentally new idea: we need to transform £1 
into a disk while controlling angular information in the Rayleigh quotients of the eigen- 
functions. Any such transformation will change the Rayleigh quotients substantially, and 
so we must devise a scheme for extracting the geometric effect and leaving behind the 
portion of the Rayleigh quotient that corresponds to the eigenfunction of the disk. 

We construct a geometric transformation that maps linearly on rays and has constant 
Jacobian. As Figure 2 indicates, wherever the transformation stretches radially it must 
compress angularly. Section 2 gives the precise definition. We will extract the geometric 
contribution to the Rayleigh quotient of the trial function on f2 by composing the trans- 
planted eigenfunction with an arbitrary orthogonal transformation U of the ball and then 
averaging over all such U (see Proposition 5.2 and Section 6). The constant Jacobian re- 
quirement is used here to guarantee that transplanting orthogonal eigenfunctions from the 
ball will yield orthogonal trial functions on VL. 

Note that Polya and Szego's result on the first eigenvalue was proved by a linear-on- 
rays transformation that does not distort angles and hence does not preserve area. In other 
words, they took 6 = <p'm Figure 2. They also did not average over rotations. Instead they 
relied on the special fact that the first Dirichlet eigenfunction of the disk is radial. 
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Perturbations of the disk. To make the theorem more concrete, we examine the case of 
nearly circular domains. Suppose P(9) is a Lipschitz continuous, 2n -periodic function. 
Define a plane domain f2 £ = {re ld : < r < 1 + eP(9)}, and assume e is small enough 
that the radius 1 + eP(9) is positive for all 9. We may regard f2 £ as a perturbation of the 
unit disk D, since f2 = D. Let j oi ~ 2.4 be the first zero of the Bessel function J , and 
recall that Ai(D) = jl x . 

Corollary 1.2 (Nearly circular domains). The first eigenvalue of the domain Q E is bounded 
above and below in terms of the boundary perturbation P: 

1< Mwi< 1+ ^r w * (2) 

JoV Jo (l+eP(9)) 2tt 



as e — )■ with P fixed. 

Eigenvalue sums satisfy a similar upper bound, for n > 1 and s £ (0,1]: 

(ELiAjW) 1 ''^) 



(E;=iA,(D) s rA(D) 

< max < 1 + 



2n jo.m2 de ^ (l+eP(9)fd9/2ii 



(l + eP(9)) 2 2^ [J^(l + £ P { 9)fd9/2n] 



1 + 0(e 2 



The lower bound on the first eigenvalue in (2) is the famous Faber-Krahn inequality. The 
upper bounds are immediate from Theorem 1 . 1 , by taking R = 1 + eP in the formula for Go 
in Proposition 10.1 and remembering that our two definitions of Go agree in 2 dimensions 
(Lemma 10.2). 

The upper bound on the first eigenvalue in (2) is equivalent to the estimate of Polya and 
Szego [38, pp. 14-15, 91-92]. The upper bound on eigenvalue sums is new. 

Let us compare the upper bound on the first eigenvalue with recent work of van den Berg 
[8, Theorem l(ii)]. He obtained an estimate of the form 

1 + C 2 (||P|| 2 ||P'|| 2 + \\P\\l)e 2 + GallP'H^ 3 , (3) 

for perturbations normalized by ||P||oo = 1- Note his formula involves HPl^H-P'lh at 
second order whereas our corollary has HP'Ul. When applied to the main example in van 
den Berg's paper, his estimate beats our estimate (2) by a factor of e 1 ^ 2 because his P 
depends on e with ||P|| 2 = O^ 1 / 4 ) and ||P'|| 2 = O^ 1 / 4 ). On the other hand, when 
applied to the "uniform" perturbation P = — 1, equality holds in our estimate (2) whereas 
(3) is not exact, due to the contribution of ||P|| 2 . 

Next let us contrast with Rayleigh's second order perturbation expansion of A x . See 
either his formal derivation [40, §210], or Polya and Szego's account [38, pp. 132-133], 
or Henry's more modern approach in all dimensions [18, p. 35]. This expansion gives no 
error bounds, and it holds only for P fixed with e tending to 0. Both our corollary and van 
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den Berg's work establish approximation bounds, and allow P to vary with e. Rayleigh's 
perturbation formula is better in one respect, though, because it gives an exact second order 
term (£ 2 -term) for Ai A. This term behaves like \n\ \P(n) | 2 , that is, like the square of the 
i/ 1 / 2 -norm of the boundary perturbation. In contrast, both our corollary and van den Berg's 
work have second order terms that are bigger, being controlled by the square of the H 1 - 
norm, ^n 2 |P(n)| 2 . 

Prior work, and the new methods. The idea of obtaining eigenvalue bounds by trans- 
forming a domain and averaging over rotations appeared already in Laugesen and Mor- 
purgo's 2-dimensional conformal mapping approach [32]. Their averaging task was much 
easier, though, since it needed only subharmonicity of the modulus of an analytic function. 
Further, their "reciprocal eigenvalue" results are inherently less powerful than the methods 
of this paper since, for example, they cannot yield the heat trace for all t > 0. 

More recently, sharp eigenvalue bounds on linear images of rotationally symmetric do- 
mains (such as regular polygons) were obtained by the authors and collaborators [30, 
31, 33, 34]. For example, they showed that the centered equilateral triangle maximizes 
(Ai + • • • + \ n )A/Gi among all triangles. The averaging in those papers takes place over 
discrete groups of rotations (such as 3-fold rotations for triangles), and relies on "tight 
frame" identities which are special cases of Schur's Lemma from representation theory. 
The transformations in those papers are globally linear, and so are simpler than the linear- 
on-rays transformations constructed in this paper. This simplicity comes at the cost of a 
severely restricted class of image domains. 

We must push beyond these existing averaging methods, in this paper. One serious 
obstacle is the nonlinear nature of our transformation, which causes the rotation matrix U 
to appear multiple times in the transformed Rayleigh quotient, both inside and outside the 
derivative of the transformation. We describe how to overcome these obstacles in Section 6. 

Complementary bounds in the reverse direction. Rayleigh conjectured in 1877, and 
Faber and Krahn proved in the the 1920s, that 

XiA is minimal for the disk. 

Many proofs and extensions are known [4, 7, 24, 38]. The result holds trivially for the first 
Neumann eigenvalue which equals zero for each domain. The result holds also for the 
first Robin eigenvalue, assuming a positive Robin parameter, by work of Bossel and Daners 
[9, 13] that was improved to irregular domains by Bucur and Giacomini [11]. 

This Rayleigh-Faber-Krahn inequality does not extend to sums and products of eigen- 
values, as the results in this paper do. It does extend to the spectral zeta function and heat 
trace, for s < — 1 and t > respectively, each of them being maximal for the disk of the 
same area under Dirichlet boundary conditions. This important extension is due to Lut- 
tinger [35], whose multiple integral rearrangement techniques proved remarkably fertile 
in joint work with Brascamp and Lieb [10]. Note that letting t — > oo in the heat trace 
inequality yields a new proof of the Faber-Krahn theorem. 

Luttinger's methods and results do not extend to Neumann boundary conditions. In fact, 
for large t the area-normalized Neumann heat trace equals approximately 1 + e~' M2At , which 
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is minimal for the disk (not maximal) by Szego-Weinberger's result on the second Neu- 
mann eigenvalue [41]. One naturally conjectures that the Neumann heat trace is minimal 
for the disk of the same area, for each t > 0. This problem remains open. For the analogous 
problem of the heat trace on the sphere, Morpurgo [37] has proved local minimality at the 
round metric. 

To put this paper in context, then, one may imagine introducing the geometric factor 
G into the classical Faber-Krahn and Luttinger inequalities in order to obtain inequalities 
in the reverse direction. The resulting theorems improve in three respects on the classical 
antecedents, because they hold: for finite sums and products of eigenvalues, for each partial 
sum of the spectral zeta function and heat trace, and for Neumann and Robin boundary 
conditions in addition to Dirichlet. 

2. The volume preserving transformation, and geometric factors 

Write B = M d for the unit ball centered at the origin, § = § d_1 for the unit sphere, and 
M + = (0, oo) for the positive half-axis. 

We say a domain in R. d is Lip schitz- starlike if it can be expressed in the form 

n = {r£ : £ G S,0 <r < R(£)} 

for some positive, Lipschitz continuous function R(£) on §. Call R the radius function of 
£1. The gauge function is its reciprocal, 

Write V for volume in R d . Notice V(fy = \ f s R(£) d dS(£). 

The volume preserving (constant Jacobian) transformation. Our work relies on a map 
from the Lip schitz- starlike domain fi to the ball that preserves volume locally (up to a 
scale factor), and is linear on each ray from the origin. See Figure 2 for an example in 2 
dimensions. First we construct a "boundary homeomorphism" H associated with the radius 
function R of fi. 

Lemma 2.1. There exists a bi-Lipschitz homeomorphism : § — > S> that distorts surface 
area in proportion to the d-th power of the radius function: 

JacMO = £gfl({)< (4) 

We prove the lemma in Section 11. Simply note at this stage that the left and right sides 
of (4) both integrate over § to yield \S\. 

Now define the mapping T : Q — > B by mapping linearly in each direction and trans- 
forming directions with H; that is, define 

T(r « = W) m) <5) 

for vectors £ G § and numbers r G [0, R(£)). One can check that T is a bi-Lipschitz 
homeomorphism of £1 to B. Its Jacobian determinant is constant, with 

Jac(T) = V(M)/V(Q), Jac(T _1 ) = V(n)/V(M), 
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as one deduces from the definition (5) and Jacobian formula (4). 

The constant Jacobian property of T will be essential later, when we transplant an or- 
thonormal collection of eigenf unctions on the ball to a collection of functions on f2. The 
transplanted functions will remain orthogonal, thanks to the constant Jacobian condition, 
and so we can use them as trial functions in the Rayleigh principle for the eigenvalue sum. 

The geometric factors. From now on, we extend R and H by homogeneity to be defined 
not just on the unit sphere but on all nonzero vectors: 

R(r£) = R(0, H{r{) = H(£), 

for all r > 0. Thus it makes sense to speak of the gradient vector Vi2, and the derivative 
matrix DH. 

Given a real matrix M, write its Hilbert-Schmidt norm as 

||M||^=(X;^) 1/2 = (trMtM) 1 / a , 
j,k 

where denotes the transposed matrix. All matrices in this paper will be real. 
Now define two geometric quantities 

r fn , _ lif/sWo^ + iv^or^- 4 ] dm 

G o^J ; ; (d _ 2)/d , (o) 



faf s R(Z) d dS{t) 



\{d-2)/d ■ (/) 



JS| 

Clearly G and G% are scale invariant, meaning that Gi(Cl) = Gi(aQ) for all a > 0, since 
aQ has radius function aR. 

Alternative formulas for G and G\ of a geometric nature will be developed in Sec- 
tion 10. There we express G in terms of a support-type integral over the boundary 
that was employed previously by Polya and Szego, and we show in two dimensions that 
G\ = 27r/ or i g in/v4 2 . Thus these alternative formulas recover the definitions (1) that we used 
in the plane, and show that in 2 dimensions, both Go and G\ depend only on the shape of 
Vt and on the choice of origin. 

In higher dimensions, G\ depends also on the choice of homeomorphism H. 

Example. If Vt is a centered ball one has R = const., so that G = 1. By convention, for 
a centered ball we choose the homeomorphism H to be the identity on the sphere, so that 
G 1 = l. 

Lemma 2.2. The geometric quantities are always at least 1 in value: 

G > 1 and G 1 > 1. 

Equality statement: Go = 1 if and only ifVt is a centered ball, and G\ — 1 if and only ifQ 
is a centered ball and H is an orthogonal transformation of the sphere. 

This lemma helps us interpret the main theorem below. We do not otherwise need the 
lemma though, and so we defer its proof to Section 10. 
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3. Main results 

First we extend the eigenvalues estimates in Theorem 1.1 to all dimensions. Let 

G = max{G , G\}. 

Theorem 3.1 (Dirichlet). Assume Q is a Lip schitz- starlike domain in IR d , d > 2. Suppose 
$ : K + — > M. is concave and increasing, and let n > 1. Then the scale invariant eigenvalue 
functional X]j=i &{XjV 2 ' d /G) is maximal when Vt is a centered ball. 
In particular, the junctionals 

XxV't/Go, X 2 V 2/d /G , (Xi + '-' + X^V^/G, 0^X 2 ■ ■ ■ X n V 2 ' d /G, 

are maximal when Q is a centered ball, for each exponent < s < 1. For s < < t the 
Junctionals 

n n 

^(A^ 2/ 7G) S and ^eM-XjV 2/d t/G) 

3=1 3=1 

are minimal when Q is a centered ball. 

Equality statement for the first eigenvalue: if X\V 2 l d /Go\ n = X\V 2 / d /Gq\ m and R is 
C 2 -smooth then Q is a centered ball. 

The proof is in Section 6. 

The only part of the theorem known previously was the estimate on the first eigenvalue. 
This extremal result for X 1 V 2 ^ d /G was proved by Polya and Szego [38, pp. 14-15, 91-92] 
in 2 dimensions, and by Freitas and Krejcifik [16, Theorem 3] in higher dimensions. The 
geometric factors in those papers look different from our G , but they are equivalent as we 
explain in Section 10. 

Note that for the first (and second) eigenvalue, the conclusion of our theorem is stronger 
than for the general case, because it uses Go instead of G. The underlying reason is that 
the first eigenfunction of a ball is purely radial, so that our proof does not depend on the 
angular information encoded in the homeomorphism H and factor G\. 

We strengthen the theorem in Section 9 by adapting the geometric factor to each eigen- 
value. There we replace G with a convex combination of G and G\ (rather than their 
maximum). 

Perturbations of the ball. Let us see what Theorem 3.1 says for nearly spherical domains. 
Suppose P(£) is a Lipschitz continuous function on the sphere §, and define a domain 

Q £ = { r £ : < r < 1 + eP(0}» assuming < 1 + eP(g) for all f . 

Corollary 3.2 (Nearly spherical domains). The first eigenvalue of the domain f2 £ can be 
bounded above and below in terms of the radial perturbation P: 

1 S ^Ta\frfTa\-?.M - ^nlK 



as e — > with P fixed, where P = J § PdS/\S\ is the mean value of the perturbation. 
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The lower bound is simply the Faber-Krahn result. The upper bound appears not to 
have been stated before. It follows by straightforward calculations from Theorem 3.1, 
simply substituting R — 1 + eP into the definition (6) of G . It can be compared with 
van den Berg's result [8, Theorem l(ii)], just like in 2 dimensions — see the remarks after 
Corollary 1.2. 

Amusingly, the corollary implies a Poincare inequality on the sphere, since the e 2 -term 
is necessarily nonnegative. 

Neumann and Robin boundary conditions. Denote the Neumann eigenvalues by p,j, 
assuming that dVl is Lipschitz so that the spectrum exists and is discrete. Write Uj for 
corresponding orthonormal eigenfunctions. Then 

{— Auj = fijUj in f2 
= on dn 

and 

= Hx < Hi < fi 3 < ■ ■ ■ ■ 
We will ignore the first eigenvalue, in the next theorem, since /xi = for every domain. 

Theorem 3.3 (Neumann). Assume Q is a Lipschitz- starlike domain in IR d , d > 2. Suppose 
$ : R+ — > R is concave and increasing, and let n > 2. Then the scale invariant eigenvalue 
functional YTj=2 ®(HjV 2 ^ d /G) is maximal when Vt is a centered ball. 
In particular, for < s < 1 the functionals 

H2V 2/d /G, (j4 + ■ ■ ■ + Hnf IS V 2/d /G, -^2-/inV 2/ 7G, 

are maximal when Q is a centered ball. For s < < t the functionals 

n n 

J2(^V 2/d /GY and J]exp(-/i i V 2/d t/G) 

3=2 3=2 

are minimal when Q is a centered ball. 

Equality statement for the first nonzero eigenvalue: if R and H are C 2 -smooth and 
f J i 2 V 2 / d /G\ n = [i2V 2 / d /G\ B then Q is a centered ball. 

Section 7 has the proof. Because G > 1, the bound on /i 2 V A2,/d / G in Theorem 3.3 follows 
from the Szego-Weinberger theorem that [i2^ 2 ^ d is maximal for the ball (see [41], or [19, 
Theorem 7.1.1]). Note that Theorem 3.3 holds for higher eigenvalue functionals too, which 
the Szego-Weinberger theorem does not. 

Theorem 3.3 might hold with G replaced by the smaller quantity G\ (the moment of 
inertia type functional), as we have conjectured elsewhere [33, §4]. 

Next we turn to Robin boundary conditions. The Robin eigenvalue problem is 

—h 2 Auj = pjUj in Q, 
h 2 ^ + auj = on dQ, 

with eigenvalues 

< p X < P2 < P3 < • • • , 
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where a G L°°(dQ) is the Robin parameter and h > is the Planck constant. Existence 
and discreteness of the Robin spectrum under these assumptions follows from the usual 
quadratic form approach; see for example [29, Chapter 5]. Note the Robin parameter a 
will be fixed in our work. For interesting asymptotic behavior of Robin eigenvalues as a 
varies and approaches ±oo, see [14, 25] and references therein. 
The Robin eigenvalues reduce to Neumann when U = 1, a = 0. 

In the Dirichlet and Neumann eigenvalue problems we took h = 1. That causes no loss 
of generality, since one can always adjust the value of h by rescaling the domain. In our 
Robin result below, though, the Planck constant and Robin parameter will be multiplied by 
different geometric factors. Accordingly we write pj = Pj(£l, h, a) to display the depen- 
dence of the jth Robin eigenvalue on the domain, Planck constant and Robin parameter. 

Our theorem will involve a new geometric factor, 



G R obin {\ dM y V ( M yd-i)/d ) ■ W 

Clearly G Rob [ n > 1 by the isoperimetric inequality, with equality if and only if SI is a ball. 
And G is larger than G Robin : 

Lemma 3.4. G > G Robin > 1. 

The lemma was proved in 2 dimensions by Aissen [1, Theorem 1]. Our proof appears in 
Section 10, and is valid in all dimensions. 

Now we can state our sharp upper bound on Robin eigenvalues. 

Theorem 3.5 (Robin). Assume Q is a Lip schitz- starlike domain in R , d > 2. Suppose 
$ : M + — >• K is concave and increasing, and let n > 1. Then 

n 

Y,HP^M 1/d /G l '\aV^/Gll in )) 

is maximal when Q is a centered ball and a is replaced by its average value. 
For the first eigenvalue one has a stronger result (with Go instead of G): 

Pl (Q, hR/Gl /2 , aR/GlUn) < Pl (B, h, a) (9) 

where R is the radius of a ball having the same volume as Q and a = j dn a dS/\dQ\ is the 
average value of the Robin parameter. If equality holds in (9) and if R is C 2 -smooth, then 
Q is a centered ball. 

See Section 8 for the proof. 

Spectral zeta function and heat trace results follow immediately, by choosing $ suitably. 
A particularly simple corollary holds for the ball: averaging the Robin parameter in- 
creases the eigenvalue functionals on a ball, with 

n n 
3=1 3=1 

whenever n > 1 and $ : R + — > K is concave and increasing. 
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Remark. Bareket [6, Appendix A] applied the Pdlya-Szego trial function technique with 
a constant parameter a < and got an upper bound on p 1 . She did not attach geometric 
factors to the Planck constant or Robin parameter, though, and so her bound looks different 
from ours in (9). 

Bareket raised an analogue of Rayleigh's Conjecture for negative Robin parameter, namely 
that 

Pi(Q,l,a) < P i(M,l,a), 
whenever a < is constant and Vt has the same volume as B. Notice the ball is the 
maximizer here, which is the opposite of the Faber-Krahn type result by Bossel and Daners 
that holds when a > 0. Let us compare Bareket's conjecture with our result (9) for the first 
eigenvalue, which says (when h = 1, R = 1 and a = a < 0) that 

Pl {Q,l/Gl /2 ,a/G^ bm ) < Pl (M,l,a). 

The factors G and G R m b are both greater than 1, and so the Planck constant and Robin 
parameter are smaller in magnitude on the left side of our inequality than on the left side 
of Bareket's conjecture. In particular, our Robin parameter is less negative than Bareket's. 
Thus while our Planck constant tends to make the Rayleigh quotient smaller than it would 
be with Bareket's Planck constant, our Robin parameter tends to make it bigger. Hence our 
result is most likely not comparable to her conjecture. 

4. Averaging over rotations and reflections 

Our proofs will involve averaging over the group of all orthogonal transformations. 
Write 7 = 7d for the Haar probability measure on the group 0(d) of orthogonal, real, 
d x d matrices. Let "Id" denote the identity matrix. 

Lemma 4.1 (Averaging in a conjugacy class). Supose M is a dx d real symmetric matrix, 
for some d > 1. Then 

[ U- 1 MUd^(U) = \tr(M)U. (10) 
Jo(d) d 

Hence for each column vector m 6 M. d , 

t(o o)(o :)^->offs). 

The lemma is a special case of Schur's Lemma from representation theory. We give a 
short proof, for the sake of completeness. 

Proof of Lemma 4.1. Denote the left side of formula (10) by 

L = f U^MUdjiU). 

Jo(d) 

For any U E 0(d), we have LU — UL by invariance of Haar measure. Let a be a real 
eigenvalue of L with eigenvector w E M. d (using here that L is symmetric, by symmetry of 
M and orthogonality of U). Then 

L(Uw) = U(Lw) = U(aw) = a(Uw). 
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Hence each vector in the orbit {Uw : U G 0(d)} is an eigenvector of L with eigenvalue a. 
The orbit spans all of W l , and so L equals a times the identity. Taking the trace yields 

ad = tr L= [ ti(U- l MU)d-f(U)= [ tr(M) drf(U) = tr(M). 

■/O(cZ) 70(d) 

Hence L equals ^ tr(M) times the identity. 

The second formula in the lemma follows immediately, by multiplying out the matrices 
and using that f Q ,„ U dj(U) = 0. □ 

5. Averaging and spherical homeomorphisms 

Consider a bi-Lipschitz homeomorphism H : § — > §, extended by homogeneity to 
R d \ {0} so that H(r£) = H(g) for all r > 0, £ G S. 

Lemma 5.1 (Orthogonality relation). For all £ G S, we /zave (DHY(£)H(£) = 0. 

Proof of Lemma 5.1. We have | 2 = 1 for all a; 7^ 0, because H takes values in the unit 
sphere. Taking the gradient of this last identity yields that 2H(xY DH(x) = 0. Applying 
the transpose and evaluating at x = £ completes the proof. □ 

The lemma implies that DH(£)y ■ H(£) = for all vectors y, so that the range of the 
derivative operator DH at £ lies orthogonal to the image vector H(£), as one would expect 
since H maps into a sphere. 

Some complicated expressions involving H can be simplified considerably, after we 
average over all orthogonal matrices U . We will encounter expressions of the following 
type when we prove our main theorem. 

Proposition 5.2. Let H : § — > § be a bi-Lipschitz homeomorphism and let £ G § be a fixed 
unit vector. 

If F G L°°(§; R d ) is a bounded, row-vector valued function on the sphere, then 
[ F(H- 1 (UC))(DHy(H- 1 (U())Ud 1 (U) = 0. 

JO(d) 

If f G L°°(S; M.) is a bounded, real-valued Junction on the sphere, then 
[ f(H-\UQ) U- 1 DH(H- 1 (UC))(DHy(H- 1 (U())Ud 1 (U) = c(Id-CC f ) 

JO{d) 

where 

Proof of Proposition 5.2. Choose an orthogonal matrix W that maps £ to the north pole, 
meaning W( = v where v — (0, . . . , 0, l)t is the north pole column vector. 
We have that 

[ F(H- 1 (UC))(DH)^(H- 1 (UC))Ud 1 (U) (11) 

JO(d) 

= [ F(H-\Uv))(DH) ] (H-\Uv))UVW ' d^(U) (12) 

JO(d) 
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by changing variable with U (->■ UVW, where V is an arbitrary matrix in 0(d) that fixes 
the north pole (that is, Vv = v). We may write 

where V E O(d-l). 

We are permitted to average expression (12) with respect to V, since formula (11) does 
not depend on V. Notice V appears only once in (12). Averaging it yields 



/ 



Vd ld . l (V)= (° \) =M/t. 



'O(d-l) 

Multiplying on the right by gives 



/ VWd ld - X (V) = vC 

JO(d-l) 



because W( = v by construction. Hence after averaging (12) with respect to V we obtain 
that expression (11) equals 

F(H- l (Uv))(DH) ] (H- l (Uu))UvC} ' d^(U) = 0, 

0{d) 

because 

{DH)\H-\Uv))Uv = Q (14) 

by Lemma 5.1 applied with £ = z/). 

Now we prove the second formula in the proposition. We find 

f f(H-\UQ) U- 1 DH(H- 1 (UO)(DHy(H- 1 (UO)Ud 1 (U) (15) 

JO(d) 

= [ f(H-\Uv))W- l (V- 1 MV)Wd 1 (U) (16) 

JO(d) 

by changing variable with U >-)■ UVW as before, where we have defined a matrix 

M = U- 1 DH(H-\Uv))(DH) ] (H- 1 (Uu))U. 

Note that Mi/ = by (14), and similarly v^M = by symmetry of M. Hence M has the 
form 

0, 

Thus using the decomposition (13) for V, we conclude from Lemma 4.1 that 



M 
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Notice formula (15) does not depend on V. Thus after averaging expression (16) with 
respect to V and using the last formula, we find that expression (15) equals 

~ w -iflda-i 0\ w wher£ ~ = _Lf f(H-\Uv))(\xM)d-y(U). 
V u u / a — l Jou) 



Since 



W' 1 f 1 ^ 1 °^jW = W-\ld-uu ] )W = Id-CC f 



we deduce that expression (15) equals c (Id — CC + ) - 
Note that 

trM = trM = tr (dH (H- l (Uv)){DH) ] (iT^C/z/))) 

= \\DH(H-\Uv))f HS . 

Hence c can be evaluated by using the equivalence between orbital and spatial averages 
(Appendix B), which gives that 



d 



■ZJ^ll^ 11 ' 1 ^) WDHiH-'iC)) \\ 2 HS dS(('). 



Lastly, changing variable with £' = H(£) shows that c equals the constant c defined in the 
Proposition. □ 

6. Dirichlet eigenvalues — proof of Theorem 3.1 

The idea is to obtain trial functions on f2 by transplanting eigenfunctions from B to 
CI with the volume-preserving map T, and then to average with respect to rotations and 
reflections of B. 

Recall that the Rayleigh quotient associated with the Dirichlet spectrum is 

f I VwP dx 
Ray [ u \ = Jn ' ' — for u G Hq(Q). 
J n u 2 dx 

The Rayleigh-Poincare Variational Principle [4, p. 98] characterizes the sum of the first n 
Dirichlet eigenvalues as: 

Ai + • • • + \ n = min { Ray[f i] + • • • + Ray[t> n ] : 

Vi,...,v n £ -f^d(^) are pairwise orthogonal in L 2 (f2)}. 

To use this principle, let ui,u 2 ,u 3 , . . . be orthonormal eigenfunctions on B correspond- 
ing to the eigenvalues Ai(B), A2(B), As(B), .... Take an orthogonal matrix U. Then define 
trial functions 

Vj = Uj o U^ 1 o T 

on the domain fi, where the transformation T : Cl — > B was defined in Section 2 . Clearly 
Vj E L 2 (Vl) since T has constant Jacobian. One can further show that Vj has weak deriva- 
tives in L 2 (f2), since Uj is smooth with derivatives in L 2 (M) and the Lipschitz continuous 
mapping T has bounded weak derivatives. Thus Vj E Further, Vj = on <9f2 

because Uj = on <9B; more precisely, Vj E Hq(£1) because Uj E Hq{ 
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The functions Vj are pairwise orthogonal, since 

/ v jVk dx = Jac(T _1 ) / UjUkdx (17) 
Jn Jn 
= 

whenever j ^ k, using here that Uj and Uk are orthogonal and T _1 has constant Jacobian. 
Thus by the Rayleigh-Poincare principle, we have 



5>(n) <£^l 

.7 = 1 .7=1 Jj 



(18) 



The denominator of this Ray leigh quotient is f n v J fix = f Q Uj(U l T(x)) 2 dx = Jac(T ^ 
by (17) with j = A;, since the eigenfunctions are normalized with f u u 2 - dx = 1. 

For the numerator of the Rayleigh quotient, we write v = Vj and u = Uj (to simplify 
notation in what follows) and express v = v(r,£) and u = u(s, () in spherical coordinates. 
Then the relation v = u o U^ 1 o T says 

v(r,Z) = u(rr{Z),U- 1 H(()), 

by recalling that T = 1/R and using the definition (5) of T. Differentiating, we find 

v r (r,Z)=r{Z)u a (rr{Z),U- 1 H{Z)), 
(V s i;)(r,£) = ru^rTiO^-'HiO^m 

where the gradients are regarded as row vectors and by Vs we mean the gradient with 
respect to the angular variables. Writing the numerator of the Rayleigh quotient in terms 
of spherical coordinates gives that 

r r fR(0 

/ \S/ V \ 2 dx= / / (vt + r-^Vsvl^r^drdSit). 
Jn Js Jo 

After changing variable with s = »T(£) G (0, 1) and using the above formulas for v r and 
Vgv , we find that 



|Vtf dx = Qi + Q 2 + Qz 



where 



Qi= f [\m 2 + \Vm\ 2 }us(s,U- 1 H(t)fs d - 1 dsR(0 d dS(0, 
Js Jo 



Q 2 = 2 / r(£)vr(0(M)W 

Js Jo 



(V s u)t( s> U^mtyu.fa U^H^s^dsR^dS^), 



Qs = I I m'lWis^-'HiOp-'DHiOfs^dsRiO'dSiO. 
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The left side of (18) is independent of U . Hence by averaging (18) with respect to 
U e 0(d) we find 



Jac^- 1 ) 



j'=i i=i 

where we must remember that "it" means Uj, in the quantities Qi, Q 2 , Q3. 

The quantity Q 1 is easiest to average because it contains only one "IP'. We have 

[ u^U-'H&fdjiU) = -j- f u s (sX) 2 dS(C) 

JO(d) PI JS 

by the equivalence of orbital and spatial means (Appendix B). Hence 

Io id ) Qi M") 



V(fl)/V{M) 
On 
ds 



J§ 



u s (s,() dS(() s d - L ds 



by the definition (6) of G . 

For Q 2 we begin by changing variable with £ = H^iJJQ. Recalling from (4) that 
Jactf(£) = fl(0 rf V(!)/V(fi), we find that 

Q2 = 2^jj\{H-\UO)VY{H-\UO){DH)^{H-\UO)U 

(v §M )t( s ,cK( s ,c) s d " 2 ^(C). 

The integrand contains U in four locations, but averaging remains feasible: Proposition 5.2 
applied with F = TVr shows that 

[ Q 2 d 1 (U)=0. (21) 

JO{d) 

For Q 3 we again change variable with £ = i/ -1 and find that 

Q 3 = ^ jf j\(H-\UC)) 2 \VMs, QU- l DH(H-\UQ) \ 2 ^ 3 ^(C). 

In this integrand £/ appears five times, after the square of the vector is expanded. Nonethe- 
less, we can average Qs with respect to U by expanding | • • • | 2 and using Proposition 5.2 
with/ = T 2 . We find that 

f Q 3dl (U) = pQ[ [ (V s u)c(ld-(?)(Vsuys d - 3 dsdS(0 
Jo(d) V W JsJo 

where 

\DH(0\\ 2 HS 



f my 



R(0 d dS(0/ I R(0 d dS(0- 
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Note that (V§u)£ = 0, because the spherical gradient V§u lies perpendicular to the unit 
vector £. Hence we deduce that 



jaccr- 1 ; 



c ( [ (V §M )(V §M ) t s d - 3 dsdS(C) 
J§ Jo 

G ^imY'L a ~ 2]Vsu ' ?dx (22) 



by definition of Gi in (7). Combining (19)-(22) now shows that 



i=i v ' j=i 



= (Y^f d J2 [(i-^ j )G (Q)+a j G 1 (Q]\JjVu j \ 2 dx (23) 

where 

f B s- 2 \V s Ui\*dx 

(The coefficient ctj G [0, 1] measures the "angular component" of the jth energy.) Next we 
estimate G and G\ from above with their maximum G, and so conclude that 



n n „ n 

Y,^(tyV{{l) 2/d /G{tt) < V(M) 2/d Y, / \V Uj \ 2 dx = Y,\ 

3=1 3=1 "' B 3=1 



2/d 



3=1 3- 

Since G(B) = 1, we have proved the theorem in the case that $(a) = a is the identity 
function. 

The theorem now follows for any concave increasing $, by Hardy-Littlewood-Polya 
majorization as in Appendix A. 

Notice that for the first eigenvalue our proof gives a stronger conclusion, namely using 
Go instead of G, because the fundamental mode ui of the ball is a radial function and so 
«i = in the argument above. Thus \iV 2 ^ d /Go is maximal for the ball. 

For the second Dirichlet eigenvalue, we may divide and multiply by the first eigenvalue 
to obtain that 

X 2 V 2 / d _ X 1 V 2 / d A 2 
Go Go Ai 

We previously showed that the first factor is maximal for the ball, and the second factor is 
maximal too, by the sharp Payne-Polya-Weinberger result of Ashbaugh and Benguria [2]. 

Particular cases. Applying the theorem with $(a) = a s , which is concave and increasing 

when < s < 1, gives maximality of (Af H h X s n ) 1/s V 2 ' d /G for the ball. The limiting 

case s I suggests we try choosing $(a) = log a, which yields maximality of the ball for 
the functional 

n 

MA^/G) = n log ( ^\ 1 ---\ n V 2 ' d /G) . 

3=1 
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When s < we can choose the concave increasing function $(o) = — a s , which leads 
to minimality of the ball for E"=iOj^ 2/a! /G9 s . And for t > we can consider $(a) = 
— e~ at , thus obtaining minimality at the ball of Y^=i exp(—XjV 2 ' d t/G). 

Dirichlet equality statement. Assume equality holds for the first eigenvalue, that is, 

\ 1 V 2 ' d /G Q \ a = \ 1 V 2 ' d /G \ B . 

By enforcing equality in our proof above, we see that the trial function v\ on Vt must attain 
equality in the Rayleigh characterization of Ai(fi), and hence must be a first eigenfunction 
for Vt. In particular this holds when the orthogonal matrix U is the identity, so that the 
function Vi(x) = Ui(T(x)) satisfies 

Aui = -Ai(fi)ui. (24) 

The fundamental Dirichlet mode U\ of the ball is radial, with U\(x) = J(\x\) for some 
positive function J, and so we have vi(x) = J(r/R(£)). That is, 

Vl (x) = J(rT(0) 

where T — 1/R. Note R is C 2 -smooth by assumption, in this part of the theorem. 
The Laplacian of v i is given in spherical coordinates by 

A Vl (x) = J"(rr(£))r(£) 2 + ^lj'( r r(£))r(£) + iAa[J(rr(0)]. 
This spherical Laplacian can be computed by the chain rule. It equals 

A s [J(rr(£))] = j"(rr(0)r 2 |vr(0l 2 + /(rr(O)rAsr(e). 

We substitute this formula into the preceding one, and make the substitution s = rT(£). 
Then the eigenfunction equation (24) reads: 

[r(0 2 + |vr(OI 2 ] J"(s) + [(d - i)v(0 2 + r(£)A s r(0]^'(«) = -Ai(n) J (a). 

Integrating over £ G § yields that 

\\t\\ 2 h1 j"(s) + [(d - i)||r|| 2 - ||vr|| 2 ] - s j'(s) = -A x (n)|s| 

The eigenfunction equation for the unit ball (the case T = 1) says that 

J"(s) + d —^J'(s) = -Ai(B) J(s). (25) 
s 

We subtract || T times this equation from the previous equation, thereby obtaining a first 
order equation for J: 

||Vr|| 2 J'(s) = -(Ai(fi)|S| - Ai(B)||r||^i) sJ(s), < s < 1. 

Suppose T ^ const., which ensures that || Vr|| 2 > 0. Then the last equation for J has 
the form 

J'(s) = asJ(s) 

for some a E K, and so J"(s) = (a 2 s 2 + a)J(s). Substituting these relations into the 
eigenfunction equation (25) implies a 2 s 2 + ad = — Ai(B) for all s £ (0, 1), and so a = 
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and then Ai(B) = 0. This contradiction tells us that T = const., and so R = const., which 
means that is a centered ball. 



7. Neumann eigenvalues — proof of Theorem 3.3 

For Neumann boundary conditions, the Rayleigh quotient and Rayleigh-Poincare Prin- 
ciple are just as for the Dirichlet case, except using trial functions in rather than 
Hq(Q). Thus we may follow the proof of Theorem 3.1, except using Neumann eigenfunc- 
tions of the ball instead of Dirichlet eigenf unctions, to prove that 

n n 
3=1 3=1 

On each side, the term with j = 1 may now be discarded because fii = 0. Then the proof 
can be completed by majorization. 

For the equality statement, rather than adapting the Dirichlet case we present a simpler 
approach. Suppose equality holds for the first nonzero eigenvalue, that is, 

toV*"/G\ a = H2V*/ d /G\ B . (26) 

Since /i2^ 2 ^| n < A i 2^ 2//d | B by the Szego-Weinberger result [41] (or see [19, Theorem 
7.1.1]), and since G(Q) > 1 = G(M), we conclude from equality holding in (26) that 
G(Q) = 1. Hence Go(Q) = 1, and so Vl is a centered ball by the equality statement in 
Lemma 2.2. 



8. Robin eigenvalues — proof of Theorem 3.5 

The Rayleigh quotient for the Robin problem is 

r n ti 2 L\Vu\ 2 dx + fr, n au 2 dS -, , 

Ray [u] = ' - 7 Jdn foru e H (fl). (27) 

J n u 2 dx 

Let {uj} be orthonormal eigenfunctions on the unit ball B that correspond to the Robin 
eigenvalues pj (B, hViM) 1 ^, aV{M) 1 ^ d ') , for j = 1, 2, 3, .... By constructing trial functions 
and using the Rayleigh-Poincare principle as in the Dirichlet case (Section 6), we find the 
following analogue of (18): 

n 

Y,P^hV l/d l^\aV^/G)l^) 

3=1 

^ h 2 V 2 ' d L \Vvj\ 2 dx " V l l d / av 2 dS 
G Lv 2 dx ^ n 1 L v 2 dx 

j=l J M 3 j=l ^ Robin J ^ 3 

where Vj = Uj o U^ 1 o T. Averaging over U G 0(d) (as explained in Section 6 leading up 
to (23)) shows that the first sum in (28) is bounded from above by 

VftV(B) 2/d / \V Uj \ 2 dx. (29) 
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For the principal eigenvalue (n = 1) this part of the argument also works with Go in place of 
G, since U\ is radial. (When finding u\ by separation of variables, the spherical harmonics 
with angular dependence cannot arise, because u\ is positive. Hence U\ is radial.) 
We will show below that averaging the second sum in (28) gives 

n „ 

YVv(B) 1/d / u)dS. (30) 
The theorem then follows, because adding (29) and (30) gives 

n 

Pi ( B > hV(M) 1/d , aV(M) 1/d ) . 

i=i 

For (30) it suffices to consider one value of j at a time, and so we consider an arbitrary 
function u G if x (ft) with L 2 -norm equal to 1, and write v = u o U^ 1 o T. The second sum 
in (28) has terms of the form 

yV* f m av 2 dS _ f m eu(U-tT(x)) 2 dS(x) 

where in the denominator we changed variable and used that T has constant Jacobian and 
u has L 2 -norm equal to 1. Averaging the right side of (31) over matrices U G 0(d) gives 
(by the equivalence of orbital and spatial means, as in Appendix B) the expression 

This last expression equals V(B) 1//d a L u 2 dS by definition of G Rob i n and a, proving (30). 

To prove inequality (9) for the first eigenvalue, apply the theorem with n = 1 and G 
instead of G (as remarked above), and replace hby h/V{M) 1 / d and replace a by a /V r (B) 1 / d '. 

For the equality statement on the first eigenvalue, one simply adapts the proof of the 
Dirichlet equality statement in Section 6. 



9. Improvement to the main results 

Our main theorems attach the geometric factor G = max{Go, G{\ to each eigenvalue. 
An inspection of the proofs yields a stronger result, in which each eigenvalue is paired 
with a smaller geometric factor arising from a convex combination of Go and G\. We state 
this improved result below, restricting for simplicity to the case of eigenvalue sums. (The 
reader can then deduce inequalities on spectral zeta functions and so on, by applying the 
majorization result from Appendix A.) To simplify the exposition we do not treat the Robin 
case. 

Define a convex combination of the geometric factors by 

G(a; ft) = (1 - a)Go(ft) + aG x (ft), a G [0, 1]. 

To choose the relevant values of a, we fix an orthonormal basis of eigenfunctions ui,u 2 ,u 3 , . . . 
of the unit ball B corresponding to the Dirichlet eigenvalues Ai(B), A 2 (B), A 3 (B), • • • . The 
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radial energy fraction of the jth Dirichlet eigenfunction is defined to be 

£ d = k( du il ds ) 2dx 

where s G [0,1] denotes the radial variable. This energy fraction can be computed explicitly 
by writing Uj in terms of Bessel functions. Obviously < ef < 1, with ef — 1 if and only 
if Uj is purely radial. 

The angular energy fraction is then 



of = !-*{>. 



For example, the principal Dirichlet mode of the ball is radial, and so = 0. 



Similarly, we may define the angular energy fraction af for the Neumann eigenfunctions 



of the unit ball. 

Theorem 9.1 (Improved inequalities). Assume Q is a Lipschitz-starlike domain in M. d . Then 
the Dirichlet and Neumann eigenvalues satisfy 

n n 

^(tt)V({l) 2/d < Xj(M)V(M) 2/d G(af ; fi), n > 1, 

3=1 3=1 
n n 

J2fi j {Q)V{Q) 2/d < ^^(B)F(B) 2/d G(af -O), n>2. 

j=2 j=2 

Proof of Theorem 9.1. See (23) in the proof of Theorem 3.1. The Neumann case is analo- 
gous. □ 

10. Properties of the geometric factors 

Gi does not depend on if in 2 dimensions. The quantity Gi defined in (7) depends only 
on R and not on H, in 2 dimensions, by the following result. 



Proposition 10.1. In dimension d = 2, 



G o = l + ^ j 2 \\ogR)'{6) 2 dd, (32) 
^ fo n R(0f dO 2nI origin 

where A is the area offl and I origin = J Q \%\ 2 dA is its polar moment of inertia about the 
origin. Further, 

L? 

G Robin = 

where L is the perimeter ofQ. 

These formulas imply immediately that G , G\, GR bi n are > 1 and are scale invariant 
with respect to dilations of the domain, in 2 dimensions. 
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Proof of Proposition 10.1. To prove the first and third formulas, simply substitute d = 2 
into the definitions (6) and (8) of G and G^ obin . 

For the second formula, when d — 2 the definition (7) of G 1 implies 

Gi= )\DH(Zt)\\ 2 HS d9 /2tt. (34) 



The homeomorphism H of the unit circle can be written H( c s f^g) = (^^J). Homo- 
geneity of if then gives F(^) = (^^j) where 6 = arg(xi + ix 2 ) = arctan(x 2 /xi). 
Calculating the derivative matrix DH( %\ ) results in 

\\DH{f^ e )f HS ^4>'{0f. 

Further, the distortion formula (4) for H says in 2 dimensions that 

Substituting the last two formulas into (34) shows that 

J 277 R(6)U6/2tt 

Now (33) follows by evaluating area and moment of inertia in polar coordinates. □ 

Expressing G as a support-type functional, in all dimensions. The geometric meaning 
of Go is highlighted by: 

Lemma 10.2 (Equivalence of our G with the definitions of Polya-Szego and Freitas-Krejciffk). 
If Q is a Lip schitz- starlike domain then 



"2tt 

G l 



Go 

Thus in 2 dimensions, 



dn x-N(x) v 'W(n)J 

G ° = ^L^m ds{x) - (35) 

Hence ifVL is convex then Go < (V/Vi n ) 2 ^ d , where Vi n is the volume of the largest open 
ball centered at the origin and contained in Q. . 

Polya and Szego's calculations already prove the lemma in dimension 2 (see [38, p. 92]), 
but the lemma is new in higher dimensions because Freitas and Krejcifik proceeded along 
somewhat different lines in their proof. 

Proof of Lemma 10.2. Our first task is to evaluate x ■ N(x) in terms of the radius function. 
The boundary of f2 is the level set {x : \x\ 2 = R(x) 2 } and so taking the gradient gives an 
outward normal vector n(x) = x — R(x)VR(x). We evaluate at x = -R(£)£ £ dVl to obtain 
n(x) = R(g)£ - R(£)VR(R{£)£). The homogeneity relation i?(r£) = R(£) implies that 
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rVi?«) = V-R(f) for each r > 0, and so = R(£)£ - VR(g). Thus for the unit 
normal N(x) = n{x)/\n(x)\ we compute 

x-N(x)- R ^ 



^R(0 2 + \VR(0\ 2 ' 

where we used that £ • Vi?(£) = (by homogeneity of i?). 

Next we need a formula for surface area element on the boundary of f2: 

ds(x) = R(0 d -WR(0 2 + \VR(0\ 2 dS(0, 

as one proves straightforwardly by parameterizing <9f2 as {x = R(£)£ ■ £ € §}. 

By substituting the preceding formulas into the formula for Go in the lemma, we see that 
it reduces to the definition of Go in Section 2. 

Now write M in for the ball of volume V- m centered at the origin, and write Ri n for its 
radius. If f2 is convex then R in < x ■ N(x) for all x E f2, as one sees by considering a 
support plane at x, and so 

1 x- N(x) 



x ■ N(x) ~ R: 



By integrating over dQ and using the formula in the lemma for G , on the left side, and the 
divergence theorem on the right side, we find that G < (V/V in ) 2 ^ d . □ 

Evaluation of G for polygons with an inscribed circle. If is a triangle with incenter 
at the origin, or more generally if f2 is any polygon with an inscribed circle centered at 
the origin, then the geometric factor Go can be evaluated in terms of area and perimeter. 
For such domains Go = L 2 /AtiA = GR bi n , as was proved by Aissen [1, Theorem 1]. To 
prove this fact observe that x ■ N(x) equals the inradius, for each point x on the boundary. 
Hence formula (35) gives that G equals L/2ti divided by the inradius, which evaluates to 
L 2 /4:nA because A = \LR in (by triangulating the domain with respect to the origin). 

Proof that G > GR bi n > 1 and Gi > 1 (for Lemmas 2.2 and 3.4). 

We have from Lemma 10.2 and Cauchy-Schwarz that 

i |ao| 2 , v(n) \(*-D/d 

" |S| J an x-N(x)dS(x) \V(Q)J K ) 

1 \dtt\ 2 /V(B)\(*-2)/d 

by the divergence theorem (noting V • x = d), and by definition of GRobin in (8). Note that 
if equality holds then x ■ N(x) is constant, by the equality conditions for Cauchy-Schwarz. 

Further, G Rob i n > 1 by the isoperimetric inequality, as remarked after (8), with equality 
if and only if is a ball. 

Hence G > G Rob i n > 1, and Go = 1 if and only if is a centered ball. 

(Aside. This inequality was established in 2 dimensions by Aissen [1, Theorem 1].) 
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Now we prove G\ > 1. By applying the quadratic-geometric mean inequality to the 
nonzero singular values of DH we deduce that 



d-i ~ V ^ vy(o) 

where the last step uses the distortion formula for H in (4). Substituting this estimate into 
the definition of Gi in (7) shows that 

by Jensen's inequality. Since J g i? d dS* = V(VL)d and |S| = V(B)d, we deduce that G\ > 1. 

If Gi = 1 then i? is constant (by the equality conditions for Jensen), and so Q is a cen- 
tered ball. Then H maps § to § with constant distortion Jac# = 1, by (4). In 2 dimensions 
that is enough to imply H is an orthogonal transformation (cf. Section 11). So suppose 
d > 3. Note the nonzero singular values of DH all have equal magnitude (and hence 
have magnitude 1), by the equality conditions for the quadratic-geometric mean inequality. 
Therefore by the singular value decomposition, DH acts as an orthogonal matrix on the 
tangent space, at almost every point of the sphere. Liouville's theorem [21] implies that 
H is a Mobius transformation that fixes the sphere. Since also the Jacobian of H equals 1 
identically on the sphere, we conclude that H is an orthogonal transformation. 



Good choices of origin. What is a good choice of origin within the domain, given that for 
Theorem 3.1 we would like to make the geometric factors Go and G\ as small as possible? 

To minimize G\ one should choose the origin at the center of mass, because Proposi- 
tion 10.1 expresses G\ in terms of moment of inertia (at least in 2 dimensions). 

The center of mass is not generally the best choice of origin for Go- For a polygon 
with an inscribed circle, Go is minimal when the origin coincides with the center of the 
circle; this observation is due to Aissen [1, §3], with the key step being the use of Cauchy- 
Schwarz as in (36). For a triangle, for example, the inscribed circle is centered where the 
angle bisectors intersect, which can be quite far from the center of mass (as happens for a 
thin acute isosceles triangle). Thus in general one cannot hope to minimize both G and G\ 
with a single choice of origin. An exception is for domains having two axes of symmetry, 
in which case both factors are minimized when the origin is at the intersection (for Gi 
because the intersection point is the centroid, and for Go by work of Aissen [1, Cor. 1,2,3]). 

The question of whether Go or Gi is larger can be subtle to resolve. For example, 
consider the ellipse with semi-axes 3 and 1, shown in Figure 3. Choosing the origin at the 
center would minimize both factors, and in fact would make them equal (as one finds by 
direct computation). Nearby choices of origin, though, could lead to either G or G\ being 
larger. Thus for domains like perturbed ellipses it is unclear which factor will dominate, 
until computations have been performed. 



25 




Figure 3 . Choices of origin in the dark shaded region give Go *^ G\. 



1 1 . Existence of homeomorphism H — Proof of Lemma 2.1 



Here we construct a bi-Lipschitz homeomorphism H : § — > § as required in Section 2. 
We apply the construction to ellipsoids in the next section, in order to better understand the 
geometric factor G\. 

Two dimensions. The Jacobian condition (4) says H'(6) = R(8) 2 A(B)/A(£l). (Here we 
regard the radius function R as being a 2ir -periodic function of the angle 6, rather than 
being a function on the circle.) To satisfy this condition, we define 



Notice H increases by 2ir each time 6 increases by 2tt, since f Q n R(u) 2 doo = 2A(Q). Also, 
H' is continuous, and is bounded above and below. Hence H defines a C 1 -diffeomorphism 
of the circle. 

Incidentally, this construction shows that H is uniquely determined on the circle, except 
for post-rotations (adding a constant to H). We need not consider reflections (changing H 
to —H) because H has positive Jacobian and so it must be orientation preserving. 

Three dimensions: the latitude-longitude construction. Let (6i,6 2 ) be the standard 
spherical coordinates, with < Q\ < n and < 9 2 < 2n. We assume H has the form 



which means that each line of latitude {6\ = const.) is mapped to another line of latitude, 
and the longitudinal position is transformed by g. We will first determine / by studying 
how the spacing between lines of latitude must be distorted, and then will determine g by 
applying the earlier 2-dimensional method on each line of latitude. 

Fix the north and south poles (meaning /(0) = 0, f(ir) = n) and require that g increase 
by 27r for each trip around a line of latitude (meaning g(-,- + 2-k) — g = 2n). Note that 



by the Jacobian condition (4). Integrating over # 2 £ [0, 2n] gives an equation involving 
only 9\\ 




H(0 1 ,e 2 ) = (f(e 1 ),g(9 1 ,9 2 )) 



(cos A) = {smf(9 1 ))f'(e l )ge 2 (e 1 ,e 2 ) 
= £ffifl(0i,0 2 ) 3 sin0i, 




(37) 




(38) 
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This equation can be solved for cos/(#i) by direct integration, using the north pole con- 
dition cos/(0) = 1. (The south pole condition cos/(7r) = —1 then follows, since we 
have J* fo*R 3 sin 6 1 d6 2 d6 1 = f s2 R 3 dS = 3V(Q).) Next, we substitute (38) into (37) 
to get an equation for g 02 that can be integrated directly to obtain g; we fix the constant of 
integration by requiring g(9\, 0) = (which means geometrically that H fixes the prime 
meridian). One checks easily from (37) and (38) that the construction gives g(9i, 2n) = 2n 
as required. 

The above construction guarantees /' > and gg 2 > away from the poles, and one can 
check that the resulting H gives a bi-Lipschitz homeomorphism of the sphere. 

Remark. The point of this section is to provide a construction of H that can be implemented 
in practical examples. Many other homeomorphisms also satisfy the Jacobian condition (4). 
See Dacorogna and Moser [12] for an account of the amazingly varied possibilities. 

Higher dimensions. In dimensions 4 and higher, one extends the 3-dimensional construc- 
tion by means of generalized spherical coordinates. Induction on the dimension provides 
the analogue of g on lower dimensional "latitudinal spheres". We leave the details to the 
reader. 

12. Ellipsoidal examples and the geometric factor G\ 

The homeomorphism H : S — > S constructed in the preceding section induces a volume- 
preserving map T : — > B, as defined in Section 2. For ellipsoids one could alternatively 
use the linear map provided by a matrix M with M(E) = B. Which of these two maps 
will give a better estimate on the eigenvalues in Theorem 3.1? That is, which will give a 
smaller value for the geometric factor Gil 

For the linear map one has 

T(rf ) = M(r£) = rM£, 



\M(\' \M(\ 
Extending R and H to be homogeneous functions gives 

Of \ M TT( \ Mx 

R[x) = , H[x) 



\Mxy v ' \Mx\ 

Somewhat tedious calculations then show that G = G x = [V(E) /V{M)} 2 l d \\M\\ 2 HS / 'd. 
Thus for ellipsoids we recover our earlier results about linear transformations [33, 34]. 
(Those papers treat more general domains than just ellipsoids, of course.) 

Let us now compare this "linear" map T with the map T constructed from "spherical 
coordinates" as in the previous section. Table 1 shows the values of G\ in the linear case, 
and also shows numerical values from the spherical coordinates construction, for various 
choices of north pole. Notice in the table that if the ellipsoid has two equal semi-axes (so 
that it is a body of revolution), and if we take the remaining axis as the north pole, then 
we obtain the same value for Gi as in the linear case. That equality no longer holds for a 
generic ellipsoid with unequal semi-axes, and in general the linear construction gives better 
results than the spherical coordinates one does. 
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semiaxes (a,b,c) 


linear: G\ = G 


North a 


North b 


North c 


(1,1,1) 


1 


1 


1 


1 


(1,1,2) 


1.19055 


1.24002 


1.24002 


1.19055 


(1,2,2) 


1.25992 


1.25992 


1.32057 


1.32057 


(1,2,3) 


1.49810 


1.51620 


1.73826 


1.53697 



Table 1 . Values of Gi for the "linear" construction of T, and for various 
choices of north pole in the "spherical coordinates" construction of T. 




Figure 4. A sloshing cylinder C with height L and starlike cross-section fi. 

These observations provide some guidance as to how to choose the north pole when 
constructing H by the spherical coordinates method, for an arbitrary starlike domain. 



13. Sloshing problem 

We finish the paper by transferring our results to "sloshing eigenvalues". On a cylinder 
C = x [— L, 0], we consider the following two eigenvalue problems: 

Au = in C, 
u = ondC\ ft, 

du 



and 



dn 

Au 

du 

dn 
du 

dn 



Xu on Q, 



in 



on dC \ ft, 
Jlu on Q. 



See Figure 4. The second eigenvalue problem describes frequencies of sloshing of a fluid 
in the special case of a cylindrical "glass" with uniform cross-sections. (See [15] for a 
historical review and [5, 20, 26, 27, 28] for recent developments.) 
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When we want to emphasize the dependence on the cylinder depth, we write the eigen- 
values as Xj(L) and Jij(L). The eigenvalues are determined by separating the vertical and 
horizontal variables. One finds 

Xj(L) = \/XjCoih.(\/XjL), Jtj(L) = y^tarm(y//ijX), 

where Xj and fij are the Dirichlet and Neumann eigenvalues of f2. The functions $d(gi) = 
y/a coth.(y/aL) and $7v(a) = y/a t&nh(y/aL) are concave increasing, and so majorization 
(Appendix A) extends our theorems on Dirichlet and Neumann eigenvalues of the Lapla- 
cian to sloshing eigenvalues. The Neumann sloshing conclusion, in 2 dimensions, is that 
the normalized eigenvalue sum 

n r~7 r~r n 

i=2 j=2 

is maximal for the disk, for each L > 0. (Notice here that the depth L -J AjG of the cylinder 
depends on the area and geometric factor of the cross-section.) One can then extend to more 
general functionals of the Jlj, by performing a second majorization. 

These methods handle only sloshing in cylindrical glasses, although one can use do- 
main monotonicity (see [5]) to obtain bounds for some other shapes of glass. It would be 
interesting to prove sharp bounds on eigenvalues of non-cylindrical sloshing regions, by 
comparing somehow with a domain having rotational symmetry about the vertical axis, as 
in [26]. 
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Appendix A. Majorization 

To extend from eigenvalue sums to sums of concave functions of eigenvalues we use: 

Proposition A.l. Assume {aj} and {bj} are increasing sequences of positive real numbers. 
Then the following statements are equivalent: 
( i) YTj=\ a j - YTj=x b j f or each n>l. 
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(ii) Yl]=i &( a j) — Y^=i ^Kfy) f or eacn n > 1 and all concave increasing functions 
$ : R + — >• R. ^ ' 

The result is due to Hardy, Littlewood and Polya [17, §3.17]. They treated decreasing 
sequences {a,} and {bj} and a convex increasing function $, which is equivalent to Propo- 
sition A.l after replacing <3>(a) with — $(— a). A comprehensive account of majorization 
methods can be found in the monograph of Marshall, Olkin and Ingram [36]. For equality 
statements in (i) and (ii), including the infinite series case n = oo, see a paper by Laugesen 
and Morpurgo [32, Proposition 10]. 

Appendix B. Orbital and spatial averages 

Equality of orbital and spatial averages on the sphere was needed several times in the 
paper. 

Lemma B.l. 

f /(i/C) *y(tf) = 4 f f(C)dS(C) 

Jo(d) PI ./;: 

for any f G L 1 (§) and each ( G §. 

Proof of Lemma B.l. The right side of the formula equals p- f s f(U (') dS((') by a change 
of variable, for each U. Integrating with respect to U gives (by Fubini) that 

4 [ f(C)dS(C) = 4/7 /(C/C')^W)- 

For each (' we change variable with U i-> f/V, where V is chosen so that V^C' = The 
lemma follows. □ 
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